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Abstract 

We propose a definition of cycle representation for Quantum Markov Semi- 
groups (qms) and Quantum Entropy Production Rate (QEPR) in terms of the 
p-adjoint. We introduce the class of circulant qms, which admit non-equilibrium 
steady states but exhibit symmetries that allow us to compute explicitly the 
QEPR, gain a deeper insight into the notion of cycle decomposition and prove 
that quantum detailed balance holds if and only if the QEPR equals zero. 
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1 Introduction 

The notion of equilibrium state of physical systems is well understood and there ex- 
ist several conditions that characterize such states, detailed balance and zero entropy 
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production among them. For classical Markov chains the equivalence of these two 
equilibrium criteria has been proved by Qian et al.[12] using Kalpazidou's cycle repre- 
sentation for Markov chainsjH]. Non-equilibrium state is a much more subtle notion, 
since there are a huge variety of behaviors involved in it. 

This work is aimed at contributing to the program outlined in Reference[T], namely, 
to look for some interesting Gorini-Kossakowski-Sudarshan and Lindblad (GKSL) gen- 
erators with properties that are rich enough to go beyond the equilibrium situation, but 
concrete enough to allow explicit study and, in some cases, explicit solutions. We de- 
fine Quantum Entropy Production Rate (QEPR) for qms in terms of the p-adjoint 
and discuss its connection with Fagnola-Rebolledo's[6] definition. We propose a defini- 
tion of cycle representation for Completely Positive (CP) maps and GKSL generators, 
discussing its connections with the QEPR. To test and illustrate the above notions, 
we introduce the class of circulant qms that admit non-equilibrium steady states but 
exhibit pretty symmetries which allow us explicit computation of the QEPR. The sym- 
metry properties of our semigroups arise from an abelian group structure on the state 
space of the associated classical Markov chain. 

Section 2 is a brief review of quantum detailed balance and its extensions. Our 
QEPR definition along with some basic properties are discussed in Section [3j A brief 
review of cycles and passage matrices is made in Section H] . In Section |5] we show 
how the transition probability matrix of a Markov chain on a finite abelian group is a 
circulant matrix, which is the leading concept of this article. Section |6] offers a quantum 
generalization of the former, named circulant operator, and we use it to define circulant 
qms; also in this section we propose a definition of quantum cycle representation for 
CP maps and GKSL generators. In Section [7] and [8], both QEPR and classical EPR are 
explicitly computed and compared for a circulant qms and its classical restriction using 
a diagonal invariant state. The remaining invariant states and its QEPR are studied in 
Section [9] . 

2 Preliminaries 

2.1 Quantum detailed balance 

For uniformly continuous qms on 13(h), with h a separable Hilbert space, a notion 
of detailed balance was introduced first by Alicki[3] and Frigerio-Gorini-Kossakowski- 
Verri[9J. Indeed, a qms with GKSL generator £ satisfies a quantum detailed balance 
condition in the sense of Ref.[3l[9] with respect to a stationary state p (i.e., tr(pC(x)) = 
0, V x G 13(h)), if there exists an operator C on 13(h) and a self-adjoint operator K on 
h such that for all x, y G B(h) the following relations hold: 

tr(pC(x)y) =tr(pxC(y)), 
C(-)-C(-) = 2i[K,-}. 
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The operator £, is called the p-adjoint of C. For a wide class of GKSL generators, 
including those deduced from the stochastic limit of quantum theory, the p-adjoint 
coincides with the time-reversed generator if quantum detailed balance holds. There- 
fore, C can be considered as an extension of the time-reversed GKSL generator to the 
non-equilibrium situation and we expect that simple non-equilibrium situations should 
appear when studying the difference between £ and £, see Accardi-Fagnola-Quezada[lJ 
and the references therein. 

Other notions of quantum detailed balance have been introduced by Fagnola and 
Umanita [H [5]. The main idea is to separate the invariant state p into two pieces or, 
equivalently, define the p-adjoint using the inner product {a,b) s = tr(p 1_s a*p s 6) for 
< s < |, and replace relations ([I]) by 

tr(p 1 - s C(x)p s y) = tr(p 1 - s xp s C(y)), 
C{.) -£(•) = 2i[K r ). 

Due to the non-commutativity, these two definitions are not equivalent in general. 
Clearly, detailed balance in the sense of ([T]) corresponds with the case s = in 

Notice that if p is a stationary state for C, hence with y = 1 in (pQ) and using that 
C(l) = we get 

= tr(pxC(t)) = tr(pC(x)), V x G B(h). 
Therefore, p is a stationary state also for L. 

2.2 The />-adjoint and special representations 

The p-adjoint (with s = 0) C of a GKSL generator £ is a GKSL generator if and only 
if the last one commutes with the modular automorphism of p, i.e., C o cr_j = cr_.j o £, 
where o\_j(a) = pap -1 , see Theorem 8 in ReferenceJS]. 

The Markov generators can be written in the standard Gorini-Kossakowski-Sudarshan 
and Lindblad (GKSL) representation 

£{x) = i[H, x] - - ^2 ( L *k L k% ~ 2L* k xL k + xL* k L k ) , (2) 
fc>i 

where H, L k 6 B(7i) with H = H* and the series J2k>i ^k^k is strongly convergent. 

Given a normal state p on B(H), a GKSL representation (j2J) of C by a bounded 
self-adjoint operator H and a finite or infinite sequence (Lk)k>i of elements of B(H) 
such that: 

(i) tr (pLfc) = for each k > 1, 

(ii) J2k>i L*kLk is a strongly convergent sum, 



3 



(iii) if X]fe>o \°k\ 2 < 00 an d c o + Z]fe>i c kL k = for complex scalars (c fc ) fc > , then c k = 
for every A; > 0, 

is called special. See Theorem 30.16 in Parthasarathy's book[10] for a proof of the 
existence of these class of representations. Special representations are unique up to 
unitary transformations. 

2.3 Weighted detailed balance 

The notion of weighted detailed balance introduced in Reference jl], was aimed at char- 
acterizing a class of GKSL generators with properties rich enough to go beyond the 
equilibrium situation but concrete enough to allow explicit study. In terms of special 
representations, weighted detailed balance is stated as follows. 

A uniformly continuous quantum Markov semigroup (7t)t>o satisfies a weighted de- 
tailed balance condition with respect to a faithful invariant state p, if its generator C 
has a special GKSL representation by means of operators H, such that here exists 
a sequence of positive weights q := and operators H',L k of a (possibly another) 
special representation of C such that the difference C p — C has the structure 

C p -C = -2t[K,-] + U, (3) 

where K = K* is bounded and 

n(x) = J2(lk-l)L':xL' k . (4) 

k 

Quantum detailed balance holds if and only if q^ — 1 for all k. 

3 Quantum Entropy Production Rate for quantum 
Markov semigroups 

In this section we introduce a notion of Quantum Entropy Production based on the 
concept of p-adjoint. As well as detailed balance, our definition depends on which p- 
adjoint is used. Our definition is slightly different from the one introduced by Fagnola 
and Rebolledo[]. Both definitions coincide in the class of circulant quantum Markov 
semigroups introduced in Section [6] below. 

Assume that C and its p-adjoint £, are GKSL generators of strongly continuous 
qms T and T, respectively, with an invariant state p. Let % t and % t denote the 
corresponding pre-dual semigroups. 

Definition 1 For every t>0, letQ t andVt t be the states (density matrices) on S(h®h) , 
with h a separable Hilbert space, given by 

n t = (i®% t )(\n p )(n p \) 
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and 

a,= (i®T«)(K){si p \): 



where 0, = 18 «.) € (ft® ft), with the ortkononnal basis of p in 

i 

h . The Quantum Entropy Production Rate of the uniformly continuous qms %, with 
respect to the invariant state p, is given by 

at t=0 

where the relative entropy of the states rj and p is defined as 

S(r], p) = tr (ri log rj — rj log p j 

if the nullspace of 77 contains the nullspace of p and 00 otherwise. 

As a consequence of Klein's Inequality, see the work of B. Ruskai [TT] . the relative 
entropy of every pair of states p, n is non-negative 

s(v,p)>o. 

Moreover, equality holds if and only if rj = p. 

In the remaining sections we compute explicitly the Quantum Entropy Production 
Rate for circulant qms. 

Remark 2 (i) In the finite dimensional case fit (resp. Vlt) is the so called Jamiolkowski^, 
or Choi-Jamiotkowski, transform of the CP map % t T p (resp. %t T p ), with 
T p (x) = p^xp^. 



(ii) A simple computation shows that try\{l p )({l p \j = tr(p) = 1, hence \Q p )(Q p \ is a 
state on £>(h £g> h) and fit is well defined. 

(Hi) In comparison to Fagnola-Rebolledo's definition of entropy production rate, we 
remark that in our definition, the Jamiolkowski transform is not modified by an 
anti-unitary operator. Moreover, instead of forward and backward two-point states 
we use as forward dynamics the time- dependent state generated by Jamiolkowski 
transform of the semigroup (%:t)t>o a> n d as a backward dynamics the one associated 
with its p-adjoint (%t)t>o- 
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4 Cycles and passage functions 



Let S be a numerable set and c a periodic function from Z into the set S. Following the 
notations of Qian et al.|12]. we call the values c(n) of c vertices (or nodes) of c, while 
the pairs (c(n),c(n + 1)) are called edges (directed edges or directed arcs) of c. The 
period of c is the smallest integer p such that c(n+p) = c(n) for all n G Z. Two periodic 
functions c and d are equivalent if one is a translation of the other, i.e., there exists 
i 6 Z such that din) = cin + i). The above is an equivalence relation and clearly two 
equivalent periodic functions have the same vertices and period. A directed circuit 
is an equivalence class of the above defined equivalence relation. Any directed circuit 
is determined either by its period p and any (p + l)-tuple (i , i l5 • • • , i p ) with i p = i ; 
or by its period p and p ordered pairs (i , ^i), (h, h), • • • , (i P -i, i P ) with i p = i , where 
%i = cin + 1 — 1), < ii < p — 1 for some n£l 

Definition 3 The cycle (or directed cycle) associated with a given directed circuit c = 
(io,ii,--- ,i p -i,i\), p > 1, with distinct vertices %q, i\, • • • i p -i, is the ordered sequence 
c = (io,H, ■ ■ - ip-ij- 

Every cycle is invariant under cyclic permutation of its vertices. We also use the notation 
c = (c(0), c(l), • • • , c(p — 1)) for the cycle associated with the directed circuit c = 
(c(0), c(l), • • • , c(p — 1), c(0)) of period p, and use the symbol c for both the directed 
circuit and the cycle when no confusion is possible. For every directed circuit c = 
(i , i\, . . . , i p -i, io), the reverse circuit c_ is defined as c_ = (i , i p -i, . . . , %\, 
When all points of c are distinct except for the extremes, then 

,^ \ f 1 if (i,j) is an edge of c; 
c{i,J) - I q otherwise. 

Hence, every cycle c has associated an unique matrix J c = (J c (i,j)), in some com- 
plex matrix space, called passage matrix of c. 



Example. If c = (0123) and c x = (0312), then 



J, 



CO 



( 1 \ 

10 

1 

\ 1 / 



and J, 



Cl 



/ 1 \ 

10 

10 

\ 1 / 



The passage matrix J CQ of the full length cycle cq = (0, 1, 2, . . . ,p — 1) is often called 
the primary permutation matrix and the cycle Co the primary cycle. From now on J p 
will denote the primary permutation matrix in a p x p complex matrix space. Notice 
that given any cycle c = (c(0), c(l), • • • , c(p — 1)), its passage matrix can be written in 
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terms of the canonical basis {|ei)(e 3 -| : < i,j < p — 1} of the p x p complex matrix 
space as 

^ \ e c(i))(e c (i+i)\, 



Jr. 



i=0 



where {ej}o<j< P -i is the canonical basis of C p . Notice that J c moves the canonical basis 
of C p according to the cycle c, i.e., J c e c (i) — ^cU-i) f° r an So, the primary permutation 
matrix J p is, in fact, the left shift operator for the canonical basis in C p . 

5 Circulant matrices 

5.1 Markov chains on finite groups 

Let (G, o) be a finite group. Unless otherwise specified, we let p — \G\ and denote by 
hg the product h o g, h, g G G. Given a probability distribution /i on G, the transition 
probabilities 

P(g,hg) = fi({h}), 
define a discrete time Markov chain on G. 

Example 1. Consider the cyclic group Z p = {0, 1, • • • ,p — 1} and any distribution 
probability a = {ao, a\, • • • , a p _i} on Z p . Then the transition probability matrix is the 
circulant matrix 



A = circ(a , cei, a 2 , • • • , a p -i) 



a P - 2 oip-i a 



a 2 a 3 



OLr, 



Op- 2 
"p-3 



«0 / 



Notice that A is a convex linear combination of powers of the primary permutation 



matrix J p = X^=o l e i)( e i+ih indeed, 



p-i 



(5) 



i=o 



Example 2. Let G be the abelian group G = Z p x Z q where the symbol x denotes 
direct product, with p, q > 2. We set the lexicographic order in Z p x 7L q and take 
a = {a(0,0),--- ,a(0,g-l),a(l,0), • • • , a(l, g- 1), • • • , a(p- 1, 0), • • • , a(p- 1, g- 1)} 
any probability distribution on G. One can easily see that the corresponding transition 
probability matrix is the block circulant matrix 



R — circ(i?o, Ri, • ■ ■ , Rp-i), 
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with circulant blocks 

/ a(i,0) 
a(i,q- 1) 
a(i,q- 2) 



a(i, 1) a(i,2) 
a(«,0) a(i, 1) 
a(i, q — 1) 0) 



V 



a(i, 3) 



a(z, g — 1) \ 
a(z,g-2) 
a(i, q — 3) 

a(«,o) y 



0, 



The above matrix R is a convex linear combination of tensor products of powers 
of the primary permutation matrices J p = Yli=o l e i)( e «+il an d J q = Sj=o l e j)( e i+il) 
indeed, 



R 



E 

0<i<p-l,0<?'<(j-l 



(6) 



Remark 4 Z>ue to Birkhoff's Theorem, every bi-stochastic matrix is a convex linear 
combination of permutation matrices. Notice that |3]) and are Birkhoff's represen- 
tations of bi-stochastic circulant and block- circulant matrices, respectively. 



5.2 Diagonalization of circulant matrices 

The discrete (or quantum) Fourier transform on C p is the unitary operator defined by 
means of ^ 

* 0<j,fc<p-l 

where u p is a primitive p-th root of unity and {ej}o<j< P -i is the canonical basis of 
C p . Before proving the discrete Fourier transform diagonalizes J p we will need the next 
ortoghonality relation between the p-th roots of unity 

Proposition 5 For every pair i,k e {0, 1, . . . ,p — 1} 

1=0 

Proof. Fix i, k e {0, 1, . . . ,p — 1}. Any primitive p-th root of unity w p satisfies 



p Pp L^i p p p L^i p p p L^i p p ' 

therefore (u p k ^ — 1)~ ^^u^'^ 1 = 0. Since ^ — 1 = 1 — 5^ the conclusion follows. 



□ 
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Lemma 6 Let Z p = diag(l,u p ,Up, ■ ■ ■ ,oJP x ) ; then 
(y FpJpFp Zp, 

(ii) (F p g> F q ){J p ® J q ){F p g> F q )* = Z V ® Z q . 
Proof. Direct computations show that 

FpJp = — yV'|e fe )(e m | 



- E (E^'H^^ 1 = E^i e *>< e *i = z *- 



This proves (i). Item (ii) follows directly from (i). □ 
Since each circulant matrix can be expressed in terms of the primary permutation 
matrix J p , it follows that the discrete Fourier transform diagonalizes every circulant 
matrix as well as block circulant matrices with circulant blocks. 

Theorem 7 If A = J2i a (^)^p an( ^ B = J2i j a (hj)Jp ® Jq, then 
(i) 

F p AF; = J2 X k\e k )(e k \, 

k 

with \ k = ^ a(i)uJ?f , and 

(ii) 

(F p ® F q )B(F p ® F q )* = J2Ki\ 

k,l 

with X H = J2ij <*(i, j)u£u? q l ■ 
Proof. Using the above Lemma El a direct computation shows that 



f p af ; = zZ a ^ F p Jl P F P = E^E^i^i = E (E a ^?) |efc)(efel 



k i 



This proves (i). 

Now observe that 
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(F p ® F q )B(F p ® = J2 a(i,j)(F p j;F;) ® (F.J^*) 

= E^^')(^^ |efe)(efc ) (E^>h 

i,j k I 

= E ( E a&fitffflf) \ek®e l )(e k ®e l \. 

k,l i,j 

This finishes the proof. 

Corollary 8 With the notations in the above theorem we have 
(i) 



j,i 

with $ m (t) = E fe ^ fe e' Afe , and 

e tB = J_ ^ $ m _ ijn _j(t)|ei <g> e.,)(e m <g> e n |, 

i,j,m,n 

with $ y (i)=E M <^e^. 
Proof. The result of the above theorem and a direct computation show that 

F;diag(e^)F p = W (E^'^) l e i>< e 'l = ;£*i-;(*)|e;><e,|. 

Lai i- r ■ j 



3,1 k j,l 



This proves (i). 

In a similar way we see that 



e tB = (F p ® F 9 )*diag(e afc 0(^ ® ^) 

= - E ^'e* Afci ^ m <|e,®e J )(e m ®e n | 



- E $m-i,n-j(i)|ei<8) ej)(e 



pq 

i,j,m,n 
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6 Circulant quantum Markov semigroups 



6.1 Circulant completely positive maps 

We first recall that apxp complex matrix A is called reducible if there exists a 
permutation matrix P such that 

™"=(o d)< 

where B and D are square matrices of order at least 1 . A matrix is called irreducible if 
it is not reducible. It is well known, see for instance Theorem 5.18 in Zhang's book|13J, 
that every irreducible pxp permutation matrix A is permutation similar to the primary 
permutation matrix J p , i.e., there exists a permutation matrix P such that A = PJ p P~ x . 

Lemma 9 For every irreducible permutation matrix J G A4 P (C) there exists a unique 
cycle c of maximal length such that J = J c is the passage matrix of c. 

Proof. Being an irreducible permutation matrix, J is permutation similar to J p , i.e., 
there exists a permutation matrix P such that J = PJpP^ 1 . Therefore for any element 
of the canonical basis {e{\ of C n we have 

JPc-i PJ p e\ Pei—\. 

Define a unique cycle c by means of the permutation P taking e c u\ = Pti-\. Clearly 
J = J c since Je cW = e c( j_i). □ 

Lemma 10 Let B , B\, . . . B p _i be p- dimensional subspaces mutually orthogonal with 
respect to the Hilbert- Schmidt inner product in Ai p (C), with B the subspace of all 
diagonal matrices. If J c is the passage matrix of any cycle c of maximal length, then 

J c B, = £, +1 ,V0<Z<p-l^ 

Bi = span{\e c ( k) )(e c ( k+l )\k = 0, 1, . . . , n - 1} V < / < p - 1, 
where the sums in the indices k,l is modulus p. 

Proof. Clearly condition on the left hand side of (I7|) is sufficient for J c Bi = Bi + % for all 
< I < p— 1. Let us proof the necessity by induction on p > 1. For p = 1 the condition 
on the left hand side clearly holds. Now, assuming that the condition holds for any 1 < p 
and let us proof that it holds for p + 1. We have Jc\e c (i)) ( e c(i+i)\ — \ e c(i-i)) ( e c(i+i)\ G 
Bi + i, W0<i<p,0<l<pbj assumption, where the sums in c{i — l),c(i + /) is 
modulus p + 1. Hence we have that Bi = span{|e c (j))(e c (j + i)| : 0<i<p},W0<l<p 
since these p+1 vectors are linearly independent and Bi is (p + l)-dimensional. □ 
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Definition 11 A linear operator $ : jM p (C) — > _M P (C) is called circulant map (or 
circulant quantum channel) if there exist p- dimensional subspaces B\, I = . . . ,p — 
1, mutually orthogonal with respect to the Hilbert- Schmidt inner product with Bq = 
span{\ek) (ek\ : k = 0, . . . ,p — 1}, invariant under the action of such that 

(i) M p (C) = ®^B l . 

(ii) there exists an irreducible permutation matrix J G Ai p (C) such that JBi = Bi + i 
sum modulus p. 

(Hi) If c is the cycle associate with J by Lemma\f^ then under the isomorphism from 
Bi into C n defined by \e c ^)) (e c ( k +i)\ e c ^) we have that 

$(|e c (fe))(e c (fc+i)|) !->■ e c{k) Q 
where Q is a circulant p x p matrix. 

Example. Let c be any cycle of maximal length in {0, 1, • • • ,p — 1}, then the CP 
linear map defined by 

*e(s) = J>(p 

k=0 

for some 7(7) > is a circulant CP map. Let us define the subspaces Bi = span{|e c (fc))(e c (fc + /) 
k = 0,1, . . . ,p — 1} for I = 0, . . . ,p — 1, clearly condition (i) in the above definition 
holds. With J = J c in the above definition, let us prove the invariance of the subspaces 
Bi's. For any = Y7j=o x jJ+i\ e c(j)) ( e <j+i)\ e B l we nave that 

p—i p—i 

®c(x (l) ) = ^2^l(p- k)x jJ+ i\e cij _ k) )(e c{j+ i_ k) \ 

k=0 j=0 

p-1 /p-1 \ 
= ^2[^2l(p- k )xj+k,j+k+i J |e c (j))(e c(j+0 | G Si. 

j=0 \fc=0 / 

Moreover, using the isomorphism induced by the cycle c we get 

p-i 

®c(\e c (j)}{e c (j+i)\) = ^2l(p~ k)\e c{j -k))(e c (j-k+i)\ 1 — ► 

k=0 







( 7(0) 


7(1) 


7(2) • 


• lip- 1) \ 


p-1 




lip- 1) 


7(0) 


7(1) " 


• lip -2) 


^2i(p- k ) e cu-k) 

k=0 


= e c(j) 


lip ~ 2) 


lip- 1) 


7(0) • 


• lip ~ 3) 






V 7(1) 


7(2) 


7(3) ■ 


' 7(0) ) 
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Theorem 12 For every circulant CP map $ on Ai p (C) there exists a cycle c of length 
p such that <3> = $ c . 

Proof. Assume that $ is a CP circulant map. By Lemma ( JTU1) and condition (Hi) we 
have that $(l e c(j))( e cO'+ol) = Yli=o P( k )\ e c(j+k)) ( e c(j+k+i) \ with some /3(£;)'s independent 
of / and positive. On the other side, if $ c (x) = ^2 k f3(p — k)J^xJ^* we have that 
^c(\e c {j))(e c {j+i)\) = YTJoPip- k)\e c{j -k))(e c u-k+i)\ = E fc P(k)\e c{j+k) )(e c{j+k+l) \, sums 
modulus p, therefore <£> = $ c on B[ for every < I < p — 1. Hence by condition (ii) we 
can conclude that $(a;) = 3> c (x) for all x G Ai p (C) and this finishes the proof. □ 
Consider the CP map on M P (C) <gi .M 9 (C) defined by 

= £a(p - i, q - j)(T p ® ® J|)*, (8) 

with a(i,j) > for all0<2<p — 1, 0<j<g — 1 and J s , s = p,q, the left shift 
operator. 

Motivated by the above discussion, maps of the class (jHJ) will be called block circulant 
CP maps. More generally, we call block circulant CP map to any CP linear combination 
of tensor products of powers of passage matrices. Restriction of block circulant CP maps 
to invariant subspaces coincide with block circulant matrices with circulant blocks. 

Theorem 13 For every (k,l) G Z p ® Z q let B kl be the subspace of M P (C) ® Ai q (C) 
defined by 

B M = span{\ei) (e i+k \ <g> \ej)(e j+ i\ : < i < p - 1, < j < q - 1}. 

Then, 

(i) the pq- dimensional subspaces B k i are mutually orthogonal with the Hilbert- Schmidt 
product, invariant under the action ofQ* given by ® H B kl = J\A p (C)®J\A g (C) 
and moreover, 

(ii) the restriction o/$* to any subspace B k \ reduces to the action the block circulant 
matrix Q = a (h j)(J p ® J° q ), through the isomorphism from B k i onto C p ® C q 
defined by [e^ g) 6j)(ei +k ® 6j + i\ H- <E> e 3 -. More precisely, 

$*(|ei )(eio+fcl ® |e JO )( e io+il) 

5^o;(p - (i - i), q - (jo - j)) e i ® e i = ( e i ® e io)Q- 
y 
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Where Q is the block circulant matrix Q = circ(Qo, Qi, • • • , Q p ~i) with circulant 
blocks 



( a(i,0) a(i, 1) a(i,2) ■■■ a(i,q — l) \ 
a(i,q — l) a(i,0) a(i, 1) ••• a(i,q — 2) 
a(i,q — 2) a(i, q — 1) a(z,0) ••• a(i,q — 3) 



, < i <p - 1. (9) 



y «(?, 1) ct(i,2) a(i,3) ■■■ a(i,0) J 

Proof. For every fixed (io,jo) € we have that 

$*(ho>( e j +fel ® l e io)( e jo+'l) 

= a(p - i, Q - i) ® j£) \e io <g> e jo )(e io+fc <g> e jo+ j| (j* <g> 

= a ( p ~ 9 ~~ ® ejo-j)(eio-i+* ® e j -j+i\ ^ 

a(p - i, q - j)(e iQ _i <g> e jo _j) = (e io <g> e jo )Q. 

This proves that every subspace Bki is invariant. They are mutually orthogonal, since 

tr(\e i+k <g> e j+ i){ei <g> e^- j j e» <g> ej){e i+k ' ® e i+i /|j = 5*^'. 
This proves the Theorem. □ 

6.2 Circulant quantum Markov semigroups 

Consider the discrete time Markov chain on the abelian group Z p x Z q associated with 
a given probability distribution a : Z p x Z g i— >■ [0,1] with ^ . a(z, j) = 1. If we 
set a(0, 0) = 0, then the corresponding bi-stochastic circulant transition probabilities 
matrix 

can be considered as the transition probability matrix of the embedded Markov chain of 
the continuous time Markov chain with infinitesimal generator (or Q- matrix) Q = II — 1, 
where 11 denotes the identity matrix in Ai p (C) <g> Ai q (C) = Ai n (C). Clearly Q is a 
block circulant matrix with circulant blocks, we shall consider the quantum extensions, 
in pre-dual representation, 



(10) 
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and 



C,{x) = $,{x) -x. (11) 

of II and Q, respectively, with x G Ai p (C) <g) Ai g (C). Clearly is a circulant CP 
map (embedded quantum Markov chain) . We call a circulant GKSL generator and 
circulant qms the semigroup generated by £*. 

Instead of the matrices J p and J q , we can choose any pair of passage matrices J Cp , J Cq 
of cycles of maximal length in Z p and Z g respectively, having G = Z p x Z g . Even more, 
any finite number n of maximal length cycles can be chosen in Z Pfe respectively, and 
follow the computation along the same lines with G = x^~qZ Pj .. Moreover, if Z p has a 
prime order, then every power of a passage matrix J c is the passage matrix of some 
cycle Cfc^c if O^/c^l (mod p). 

Having this in mind and Kalpaziduo's cycle representation [8] of an irreducible 
Markov chain with uniform stationary measure tv = {-} and circulant generator Q, 

-Q=Y] w c J c - -I, 

p £L p 

we can regard equations ( fTOl) and ( fill as a quantum cycle representation of the circulant 
GKSL generator with cycle weights (a(i, j))(i,j)ez p xz t - This motivates the following. 

Definition 14 Given a bounded GKSL generator of the form (Ej) with a discrete spec- 
trum Hamiltonian, we call cycle representation of its embedded quantum Markov chain 

$oo = J2 L i xL ^ 

k 

to a GKSL representation of $ of the form 

i 

where for each I, ai > and Ui is a passage matrix. 

Clearly, any cycle decomposition of the embedded chain $ induces a cycle representation 
of£. 

Remark 15 Tensor product like J l p ® J 3 q are irreducible matrices. Hence by Lemma 
[21 they are passage matrices of a cycle. This shows that our definition includes repre- 
sentations of the form l[W\) and its extensions involving any finite number of cycles or 
higher order tensor products. 
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By Theorem [13] each subspace B^i is invariant under $*, and, consequently, 
also under the action of the semigroups % = (%t) t>0 generated by £*. The state 
p = i(l <S> 1) is clearly invariant for % since 

£*(p)= Yl a(p-i,q-])l® 1-1(8)1 = 0, 

(i,i)^(0,0) 

because .w/ ^ a(p — i,q — j) = l.The p-adjoint £ has the GKSL representation 

with = L*j and Ly = a(p — i,q — j)^(Jp <S> Hence = a{i,j)^(J l p ® J 3 q ). 
One can write in direct representation 

(m)#(o,o) 

Hence the difference between £ and its p-adjoint (reverse) operator looks like 

c(x)-c(x)= Hi,j)-a( P -i,q-j))(j;®Ji I yx(j;®jj) 

(M)7H0,0) 

(i,i)^(o,o) 

with q(i,j) = a(i,j)a(p — i,q — j) -1 and the Ly's as above. Therefore, the semi- 
group 7^ satisfies a weighted detailed balance condition in the sense of Accardi-Fagnola- 
Quezadapp with weights q = [q^ = a(i, j)a(p — i, q — , see equation ([3]) above. 
Consequently, by Corollary 2 in Ref . pQ , detailed balance holds if and only if 

«(p -i,Q-j) = <*(i,j), V (0, 0) ^ (i, j) e 7L V x Z g . (12) 

7 Quantum Entropy Production Rate for circulant 
qms 

Let us compute the Quantum Entropy Production Rate (QEPR) for the circulant semi- 
group % in the previous section. For simplicity we consider first the invariant state 
p — ^1, other invariant states are studied in Section [9] We know that every subspace 
Bki of Ai p £g> Ai q is invariant under the action of the elements of %. This implies that 
the states fl t and Clt are diagonal with respect to the canonical basis. 
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Lemma 16 With the notations in Section^ and Subsection^ the following hold: 

(i) for every {i',j') G Z p x Z g; using the isomorphism induced by lemma\B, we 
have 

%t(\ei <g) ej){ev ® e r \) h-> (ei <g> 6j)e tQ = — ^ $ w _ i|W _ i (^)(e m ® e n ) 

1 m ' n (13) 

^ / ^m.n.^) \ e m+i ® &n+j) ( e m+i' ® e n+i'|- 



m,n 



where 

k,l i,j 

We recall that in this case a(0, 0) = —1 and j)^(oo) = 1- Moreover, the 

functions $ m ,„(£) are real-valued, since Q and hence e tQ are real matrices. 

(ii) 

' ^ (14) 



= — y^$m,n(^)|Mmn)( 

pq 

m,n 

where u mn = ^ ^(e, <g> e 3 -) <g> (e m+i <g) e n+j; 



Proof. Item (i) is an immediate consequence of Theorem [7] and Corollary [SJ Now a 
direct computation using (i) shows that 



/ \ / > Prs&r ® e s ® C?rt+r ® Gn+s 
m,n i,j r,s 

2pg ^ — ^ 



(15) 



This finishes the proof. □ 
The subspaces are invariant also for the reverse semigroup % t . Moreover, similar 
computations yield the following. 

Lemma 17 For the p-adjoint (reverse) semigroup we have: 
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(i) 

%t{\ei® ej)(ei> <g> e f \) H> (e* <g) ej)e tQ * = -^$ m _ !in _ i (t)(e m <g> e n ) 

! . (16) 

^ — ^ $ m ,„(t)|e m+i (g> e n+j )(e m+i / (8> e n+J -/|, 



i n. n 



where Q* is the transpose of Q and $ OT)Tl = <&,p_ mig _ n . 

^ ^ ^*p- m.q- nify .^mri) ('-'m 



(17) 



i a. 1 1 



Theorem 18 Let £* fre a circulant GKSL generator of the form (TO]), £/ien £/ie Quan- 
tum Entropy Production Rate of the corresponding qms is given by 

,_ . 11 ^-^ i t \ i \\^ a(m,n) 

e p (7*,p) = -— > (a(m,n) -a(p-m, q-n)) log— -. 

2pg, t^ v „ a(p-m,q-n) 

Proof. From the above lemmata it follows that the relative entropy has the explicit 
expression, 



S(fl t , Q t ) = — ®m,n{t)log- { 



^m,n if) 



pq ' $ p _ mj(? _ n (t) 

= ^ E (*W(*) - Vm. g -n(*))^ 

For the Quantum Entropy Production Rate we have, 
e P [%,p) = lim 



t-s>0+ t 

-~y( lim ^."W- Vm, g -n(t) N ^ ^ 



But a simple computation shows that for every m, n, 

,. $m,n(*) .. // ^ e* Q -J \ (e ® e ,e m ® e n ) 
hm = hm ( e ® e , e m ® e n ) H 

= (e <8> e , Q(e m <g> e n )) = a(m, n). 
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Therefore 



. 11 v-^ / / \ / \\ , a(m,n) 

e P {%,p) = 2~A^ ( a ( m ' n ) - a{p-m,q-n)Jlog 



2pqt-^\ ' J a(p — m,q — n) 

m,n v ' 

This finishes the proof. □ 

8 Comparison to Classical Entropy Production Rate 

The Quantum Entropy Production Rate (Op) aims at generalizing the classical one, 
hence it is natural to expect that some relation can be found between them. In this 
section we compute explicitly the (classical) Entropy Production Rate for the restriction 
of Circulant Quantum Markov Semigroups to the diagonal commutative sub-algebra, 
namely B 00 , and show it actually coincides with its quantum counterpart. 

According to Qian et al. [12] , the Classical Entropy Production Rate of an irreducible 
Markov chain with intensity matrix Q = (qi^ijes and stationary measure % = (7r.j) ig s, 
over a finite state space S is given by 



e 



v 



- ijes jQji 



(19) 



By Theorem [131 the restriction of to B 0Q reduces to the action of the block 
circulant matrix Q = circ(Qo, Qi, . . . , with circulant blocks of the form (jH]) and 

a(0,0) = —1. In terms of the distribution a, each matrix element of Q is given by 
qij = «((/ — k) p , (rj — ri) q ) where for every pair < i,j < pq — 1 we write i = qk 
j = ql + rj, < k, I < p — 1, < Ti, Tj < q — 1, and for every — (s — 1) < x < s — 1 
s = p,q, we define 



, , J x if x > 
^ X)s = \ s + x if x < 0. 

Clearly the relation ( ), holds true. 

Corollary 19 T/ie Quantum Entropy Production Rate of a Circulant qms equals the 
Classical Entropy Production Rate of its diagonal-restricted Markov chain, i.e., 

i'p (7* • p) Cp. 

Proof. An application of the above formula (Tl9l) . re-ordering the sum according with 
the order of blocks and the change of variables m = (I — k) p , n = (r 3 - — rj) q yields, 
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1 1 




6p 2pq 



k,l=0 ri,rj=0 



log 




a((l - k) p , {rj - n)g) 
a (( k ~ 0p> ( r i - r j)q) 



e p (%,p). 



This proves the corollary. 



□ 



9 QEPR with respect to other invariant states 

To close the paper, in this section we compute the QEPR in any invariant state of the 
semigroup %. 

Proposition 20 Every invariant state of £* has the form 



where pij are complex numbers constrained by the positiveness of p. 

Proof. We decompose p into its mutually orthogonal components in the subspaces 
B kl , namely p = Y^uhi- Clearly £*(p) = if and only if £(p w ) = for every 
(k,l) G Z p x Z r As a consequence of Theorem [131 using the isomorphism defined 
there, each of the above conditions becomes a linear system of equations of the form 
PmQ — 0, where Q is the same circulant matrix for all systems. Any solution to these 
systems is a multiple of the identity vector, which yields the solution f )20|) . Although 
every choice of complex constants pki give a solution of £*(/?) = 0, not all of them 
give back a state p. In fact, p o = ~ so that trp = 1 while the remaining p' kl s are 
constrained by the positiveness of p. Conversely, if p has the form f)20p then £*(p) = 



Wo / 
By Lemma [6], any invariant state p can be diagonalized by the discrete Fourier 
Transform, indeed, 




(20) 




□ 




ik 
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where p ik = h ^ Piji^u^ 1 and e z = F*e h e k = F*e k . 

PI ij 

In the next computations it is understood that sums over the first coordinate of the 
tensor product go from to p — 1 and sums over the second coordinate go from to 
q — 1. We use the results and notations in Lemma [T6l 

Let us compute the state associated with % using the basis of p, {e^ g) ej}(ij)ez p xz g , 



fit = \ 6i ® e *')( e i ® ® 7^t(p2 |gj <g> ej/)(ej (g) ey|p2) 

ii'jj' 

= — p2, / oJ.,a^ l aJ^ l 'a;^a;^- J '|e„ g> e n /)(e r g> e r /| ® 7Zt(|ei (g> ^(e,- <g> e f \) 



pq irr., 



(pq) 



L_ ^ fl^ ^(n+JV)i-(n'+7V')i' (r'+K')i' x 

„\2 2^ Pii'Pjj' U p U q U P U Q X 



ii'jj' 
nn'rr' 
NN'RR' 



\e n <g> e„')(e r <g> e r /| <g> 7^ (|ejv <g> e N >)(e R (g> e/j'|) 



=— V <f> M Mt)P(7i,N,ri,N') 



x 



nn'rr' 
NN'RR' 



MM' 

fi(r, R, r', R')\e n (g> e„/)(e r (8> e r >\ <g> \e N+M ® e NI+M ')(e R+M <g> e R > +M ' 

/ ^m,m'\tj\^mm'/\U'rnm'\j 
pq ^ 

mm 

where u mm > = -j= YLwlv ^( Z > ^> L ')( e i ® e «') ® ( e £+™ ® ez'+m'), and A ^ L 

-J- V fl* Tj(l+L)iTj(l'+L')i' 



Direct computations show that the p-adjoint semigroup, with respect to any p of 
the form (1201) , coincide with 71 given by Lemma [T71 In a similar way we get 



fit / , ^t)-m,n-m'(^)|"mra')(''r 

pq ^— ' 

m,m' 



(21) 



It follows that the Quantum Entropy Production Rate in any invariant state p of the 
form ( |20l) coincides with the one given by Theorem [181 

Theorem 21 Let % a circulant qms with GKSL generator of the form fill) , then 
the following are equivalent: 

(i) % satisfies a quantum detailed balance condition with respect to any invariant 
state p of the form K2ty) , 
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(ii) a(m, n) = a(p — m,q — n) for all (m, n) G Z p x 7L q , 

(Hi) the Quantum Entropy Production Rate of % with respect to any stationary state 
p of the form equals zero, i.e., e p (%,p) = 0. 

Proof. The equivalence of (i) and (ii) follows from Corollary 2 in Ref.pQ, see (fT2|) . 
And the equivalence of (ii) with (Hi) follows from Theorem [181 D 

Remark 22 (i) TheoremsU~B and{21\ have a direct generalization to the case of any 
finite number of cycles (or cyclic factors in the abelian group G). 

(ii) We remark that in the case of a separable probability distribution a(i,j) = a p (i)a q (j) , 
a direct computation using Lemmata XTb\ 11 shows that the states Qt, &t o,re sep- 
arable. Indeed, Q t = Q P (t) <8> Q q (t), with 

M s (t) = - J2^s(i,t)\u s (i))(u s {i)\, 

i 

where & s (i,t) = V} w y e' Aj(as ), Xj(a s ) = a s (i)uf j , and u s (k) = 
j i 

~7=/] \ e n)(e n +k\, s = p,q. 
V s ^ 
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